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Abstract: We use a diffusive model for the propagation of Galactic cosmic rays to estimate the charged 
pion production in interactions with protons of the interstellar medium. Cosmic ray nuclei from proton 
to iron are considered and the corresponding contribution to the neutrino secondary flux produced as a 
result of spallation is also estimated. 
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Introduction 

We used the numerical Galactic Cosmic Ray 
propagation code GALPROP 1 [1] to reproduce 
the diffusive galactic spectrum from SuperNova 
Remnants (SNRs) and estimate the charged pion 
production in interactions with protons of the 
interstellar medium. From this estimation we 
calculate the flux of neutrinos coming from the 
decay of charged pions. 



Diffusion Galactic model 

The diffusive model assumes cylindrical symmetry 
in the Galaxy, with coordinates R and z equal to 
the Galactocentric radius and the distance from the 
Galactic plane. The propagation region is bounded 
by R = Bj, = 30 kpc and z = Zh = 4 kpc, be- 
yond which free escape is assumed. The diffusion 
coefficient is taken as (3Dq(p/ Pd) 5 , where p is the 
particle rigidity, .Do is the diffusion coefficient at a 
reference rigidity po and S = 0.6. 
The distribution of cosmic rays sources used repro- 
duces the cosmic-ray distribution determined by 
the analysis of EGRET gamma-ray data [1]. 

Nuclei with Z < 26 are injected with a rigidity 
power law spectrum of index a = —2.05, indepen- 
dently of energy, with isotopic abundances given 
by space measurements of the cosmic ray abun- 



dances [2]. The neutrino flux is calculated as the 
product of the decay of charged pions which, in 
turn, are created in collisions of cosmic-ray parti- 
cles with interstellar gas. 

The used interstellar molecular, atomic and ion- 
ized, (H2 , HI, HII) hydrogen distribution are de- 
rived from radio HI and CO surveys in 9 Galac- 
tocentric rings and from information on the ion- 
ized component. The distribution of molecular 
hydrogen is derived indirectly from CO radio- 
emission and the assumption that the conversion 
factor H2/CO is the same for the whole Galaxy 
[3]. The atomic hydrogen (HI) distribution is taken 
from [4], with a z-dependence calculated using two 
approximation at different galactocentric distances 
R [5, 6] The ionized component HII is calculated 
using a cylindrically symmetric model [7] . 

Pion production 

Pion production in pp-collisions is calculated fol- 
lowing a method developed by Dermer, which 
combines isobaric and scaling models of the re- 
action [8]. The two models work well at low 
and high energy respectively. In the energy range 
3 GeV — 7 GeV an interpolation of the two models 
is used. 

In the isobaric model the distribution of pions is 



1. http : // galprop.stanford.edu/web-galpr op/ 
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calculated by the integration over the isobar mass 
spectrum (Mi = m p + m n9 M 2 = \fs — m p ): 

WJ? [ M2 i U(E^E p ;m A ) 

r 

where and E p are the pion and proton energy 
in the laboratory system (LS), m° A is the average 
mass of the A-isobar, T is the width of the Breight- 
Wigner distribution and yfs is the CMS energy. In 
the model it is assumed that the produced A-isobar 
of mass m A has either the same direction (+) or 
the opposite direction (-) of the colliding proton in 
the CMS. The produced isobar decays isotropically 
producing a pion with the distribution: 



fir (Eft 5 Ep 5 777* A ) = 



4777^7;/^ 



(2) 



where H [x; a, b] = 1 if a < x < b 
and H[x; a, b] = otherwise, with 

a± =7^7;(l-^4) and ^ = 7a7;(1+/?a#t) 
The Lorentz factors of the forward (+) and back- 
ward (-) moving isobars are 7^ = 7 C 7^ (l±/3 c /3^) 
where 7 C = ^/2m p is the Lorentz factor of the 
CMS in the LS and 7^ = (s + ra^ -m p )/2^ r sm A 
is the Lorentz factor of the isobar in the CMS. 
The pion Lorentz factor in the rest frame of the 
A-isobar is 7^ = (m 2 A + ml - m 2 )/2m A m 7T . 

The scaling model gives the lorentz invariant cross 
section for pion production as: 
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(3) 



where 



G v ±(E p ) = (l + 4m 2 p /s)- R , 



(4) 



G n o(E p ) = (1 + 23^- 2 - b )(l - 4m^/ S ) K , (5) 
= (Ci - C 2Pl _ + C 3 pi)/^l + 4m2/ s , (6) 

^tt = + (4/s)(p5_ +772^), (7) 



^* = 2m 7r ^sj c j 7r (p 7T cosO - p c ) ^ 
" [(s — ml — m 2 x ) 2 — Amlm 2 ^] 1 / 2 ' 

is the pion polar angle in LS, A, /3, 61,2,3, R 
are positive constants and mx is the X channel of 
the reaction (pp — > 7r ± + X). 
The energy distribution of pions can be obtained 
integrating over the polar angle 



F 7r (E 7r , E p ) 



where 
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The distribution of muon neutrinos produced di- 
rectly by the decay of pions (or kaons) produced in 
a pp-collision is given by 



F(E V ,E P ) = / dEME^Ep)— = (11) 

-f 

J El 



lib U/ JU 

71 BR 
dEir F 7T (E 7TJ E P ) — 

tin (1 — 777f # 



(l-m^/M^V 



with £™ n = - r) + M*(l - r)/(4E u ), 

E™ ax = (s - M\ + M%)/(2y/s), Bii is the 
branching ratio for meson decay in muons and 

r = ml/Ml 

Muon decay 

Since the muons originated from pions are pro- 
duced fully polarized, the energy distribution of the 
neutrinos/antineutrinos in the muon rest frame is 
given by 

^^W^W^l (12) 

where x = E' v / m^ with E' v the neutrino energy 
in the muon rest frame and the polar angle be- 
tween the neutrino and the muon spin. The func- 
tions fo(x) and fi(x) are given in table 1. 

Integrating over the polar angle and trasforming 
the distribution to the LS, the energy distribution 
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Table 1 : Functions for neutrinos from muon decay 

of a neutrino in the LS becomes: 
dn 1 



dy 0,, 



My,P,*)-PoiMy>P*)] ( 13 ) 



where y = E v jE^ is the muon energy and 
momentu in the LS and Pol^ is the muon polariza- 
tion. The functions g (y, g x (y, (3^) are given 
by 



go(y,Pn) 
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(14) 



2y/x-l dx 



with Xn 



min[l,2y/(l-/3 M )]. 
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Neutrino production 

The distributions of muons from the decay of 
mesons (pions or kaons) is given by 



dn dn 



BR 



dE v dE^ (1 - ml/M^Pj 



(16) 



M 



where M is the pion/kaon mass, Pm is the 

pion/kaon momentum in the LS and BR is the 

branching ratio of the decay. 

The energy distributions of neutrinos from muons 

through the pion (or kaon) decay in the LS is given 

by: 



F(E U ,E 7T ) = / 

J E mi 

-s 
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where y E v /E^ y min = E u /E™ ax , y max = 
E u /E™ in and 

E™ n = max[m^{El-^P;)l (18) 
E™ ax = max[m^{El + (3 n P;)l (19) 



with 7^ the pion Lorentz factor in the LS and 
£"*, p* the muon energy and momentum in the 
pion rest frame: 
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Taking into account pion production through pp- 
collisions we have: 
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Fig.l shows the neutrino production cross section 
from the decay of negative and positive pions. For 
a given proton energy, we calculated the cross sec- 
tion for the muon neutrinos produced directly from 
the pion decay and for the electron and muon neu- 
trinos produced from the muon decay. 

Neutrino flux 

Using the diffusive galactic model described in 
§2, we calculated the diffusive galactic spectrum 
from SNRs. From the proton and helium spectra 
and from the ISM gas distributions, we estimated 
the charged pion production. Using the calculated 
cross section for neutrino production, we calcu- 
lated the flux of neutrinos at Earth coming from 
the galaxy bulge, for neutrino energy in the range 
1 MeV - 100 GeV. In Fig.2, the total neutrino 
flux and the contribution of each kind of neutrino 
are shown. In Fig. 3, the skymap of the total neu- 
trino energy flux in Galactic latitude and longitude 
coordinates is shown. 
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Figure 1: Calculated cross section for neutrinos 
from negative (a) and positive (b) pions decay pro- 
duced in pp-collisions. The cross section has been 
calculated for fixed value of the cosmic ray proton 
energy. 
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Figure 2: Flux of neutrinos coming from the 
Galactic bulge produced by pion decay through the 
interaction of cosmic rays with the ISM gas. The 
curves indicates the total neutrino/antineutrino flux 
(y + 17), the electron (y e + V e ) and muon (y^ + z7 M ) 
contributions and the different muon neutrino com- 
ponents, neutrinos coming directly from the pion 



decay (y^ ■ 



7 M 



muon decay {y\fi 



) and neutrinos coming from the 
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